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EVOLUTION OF PRE-MAIN-SEQUENCE BINARY STAR SYSTEMS 
By 
Jon Kenneth West 
June 1979 

Chairman: Kwan-Yu Chen 
Major Department: Astronomy 

This paper develops a quasistatic polytropic model. 
It represents the pre-main-sequence contraction phase. The 
model was applied to both components of close binary star 
systems. The initial conditions for the models were 
determined from the Roche Lobe. The period and masses for 
the systems determine the maximum stable radii for both 
stars. Both stars were allowed to contract from the ini- 
tial models. In the early stages, it may be necessary for 
a secondary to evolve along the Hayshi Track. The time 
for each contraction had to be calculated by an iterative 
procedure. The time steps were corrected for nuclear en- 
ergy sources between iterations when needed. 

Models were calculated for eleven observed semi- 
detached close binaries published in the literature. The 
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combined mass of each system's components is between 2.5 
and 6 solar masses. First, the results of the model were 
compared with the observed systems, then a determination 
was made concerning their pre-main-sequence nature and 
age. The names and ages of the three that are likely to 
be pre-main-sequence are: U Cephei (6.1 X 102 years), 

X Trianguli (4.6 X 106 years), and W Ursae Minoris (1.87 
xX 106 years). The following four systems were found not 
to fit the pre-main-sequence model in this investigation: 
TV Cassiopeia, AI Draconis, 8 Persei and TX Ursae Majoris. 
Finally, the next four systems could be pre-main-sequence 
with modifications to the model: IM Aurigae (3.5 xX 105 
years), 6 Librae (8.2 X 105 years), U Sagittae (3.78 X 


109 years), and V505 Sagittarii (1.89 X 106 years). 
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CHAPTER I 


INTRODUCTION 


The problem of the formation of a close binary star 
system has not been solved. No detailed understanding 
exists for the physical process involved. The processes 
for the formation of a single protostar are still in- 
complete. At least three factors make the formation of 
a selfgravitating protostar very difficult: conservation 
of angular momentum, magnetic fields, and interstellar gas 
dispersion. 

Conventional studies of binary evolution usually do 
not consider any mechanism of formation. Usually model 
calculations, e. g. Pacznski (1971), are started with both 
stars on the main sequence. The star that is originally 
the more massive is called the primary. In most cases 
only the evolution of this component is followed in detail. 
Mullen (1974), however, followed the detailed evolution of 
both stars in the case of a W Ursae Majoris type binary 
model. Nevertheless initial conditions were that both 
stars were on the main sequence. 

The evolution of binary stars from the main sequence 
is idedereail to that of single stars until the more mas- 


sive star expands to fill its Roche Lobe. At this point, 
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Crawford (1955) proposed that a large fraction of the mass 
is transferred to the secondary. The roles of the primary 
and secondary may actually reverse. Sahada (1962) asked 
whether such a star, originally a certain spectral tyne, 
after acquiring mass from the other star could become a 
more massive object that nevertheless would show a "normal" 
spectrum of an earlier type. 

The mass transfer, during the reversal, goes through 
a rapid thermal phase. Then the star has a tendency to 
expand slowly on a nuclear time scale because of hydrogen 
burning in the core. The nuclear transfer ohase is the 
most likely to be observed. 

Assuming that a binary can be formed, Thomas (1977) 
has stated that the prototype of close binary evolution 
with mass transfer has confronted theoreticians with an 
increasing number of problems. He considers the vost- 
main— sequence evolution to be understood in principle. 
The problems arise as soon as we investigate individual 
systems. 

Recently, three investigations indicate the need to 
study individual systems of the pre-main-sequence type. 
Some observed systems are shown to be pre-main-sequence 
(Field; 1969), some observed systems violate the require- 
ment that a zero-age contact system cannot exist (Leung 


and Wilson; 1976), and some show mass transfer rates that 
are on a thermal time scale (Hall and Neff; 1976). These 


cases are reviewed below. 
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Field (1969) found that some close binary systems 
appeared to be pre-main-sequence. Field used data for 
the contraction of single pre-main-sequence stars published 
by Iben (1965) to find the age of the stars. Both com- 
ponents of the binary were assumed to be formed simulta- 
neously as single stars on their respective Hayashi Tracks 
(see Chapter II). Four systems were determined to be pre- 
Main-sequence. They are KO Aquilae, TV Cassiopeia, WW 
Cygni, and Z Herculis. Their total orbital angular momen- 
tum are Similar to that predicted by Roxburgh (1967) for 
stellar fission. 

Leung and Wilson (1976) have presented evidence that 
zero-age contact binary systems exist. Results of their 
photometric investigation show both components of V701 
Scorpii and those of V1010 Ophinchi are in contact. Their 
locations on the Hertzsprung-Russell diagram and their 
radii both suggest that the systems are zero-age. However, 
this result presents a problem. This is in contradiction 
to the post-main-sequence model for the origin of contact 
systems presented by Lucy (1968). This model requires the 
components of a zero-age contact system with a common en- 
velope should have equal masses. However, many detached 
binaries witn a mass ratio of unity are not observed. Con- 
tact binaries have been observed to be on the zero-age 
main-sequence, but none have a mass ratio of unit. [In con- 
nection with this, Whelan (1970) showed the possibility 


that a mass ratio of unity was unstable for pre-main-sequence 


binaries. An alternative possibility is that the adiabatic 
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constants of the stellar envelopes are not equal. The later 
possibility implies short lifetimes for the individual sys- 
tems. Bierman and Thomas (1972) allowed for different adia- 
batic constants by requiring a luminosity exchange to main- 
tain contact of the system. Their model did not solve the 
problem of the light curves for W Ursae Majoris systems, 
nor did it propose a mechanism for the required amount of 
luminosity transfer. Therefore, their origin remains a 
theoretical problem. 

Hall and Neff (1976) evaluated the veriod changes of 
24 observed semidetached binary systems. Their model as- 
sumes mass transfer from the cooler to the primary compo- 
nent. The mass loss results in a period change because of 
conservation of angular momentum. Using the observed 
period changes they calculated the mass transfer rate for 
each of the 24 systems. All systems showed mass transfer 
on a thermal time scale. This time scale differs from that 
of the conventional post-main-sequence models for the forma- 
tion of a contact or a semidetached binary system. For in- 
stance, Mullen (1974) found that each component of a vost- 
main-sequence binary expands away from the main-sequence 
on a nuclear time scale. Thus, the problem of the pre- 
or post-main-sequence nature of a binary has not been solved. 
Each system must be evaluated on an individual basis. 

The problem that the writer wishes to investigate is 
the evaluation of individual systems with a pre-main-sequence 


model. The main purpose of this investigation is to find 
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possible pre-main-sequence binary systems. This will be 
accomplished by comparing the positions of the observed 
binary on the Hertzsprung-Russell diagram with the evolu- 
tionary tracks generated by the pre-main-sequence model. 
If both components appear to have the same age and fall 
near the tracks of the model, then the system is likely 
to be pre-main-sequence. The investigation will include 


the evaluation of eleven binaries each with a total mass 


between 2.5 and 6 solar masses. In all cases the secondary 
star is over-luminous. These systems are classified as 
semidetached. 


The model being considered is independent of the ori- 
gin of a binary star system. The origin may be either con- 
densation of both components as described by Wood (1962) 
or fission of a rapidly rotating pre-main-sequence core 
as described by Whelan (1970). Yamasaki (1971) studied 
the case immediately after fission where mass transfer may 
occur. As the stars contract, they must eventually detach 
themselves from their Roche Lobes. Yamasaki ends his in- 
vestigation when mass transfer stops or a single star is 
formed because of excessive mass loss. The writer assumes 
the first case as the initial configuration of the pre- 
main-sequence model. At this stage the individual stars 
fill their Roche Lobes. Each star is allowed to contract 
gravitationally and independently toward. the main-sequence. 

As a comment on the study of origins, I would like to 


quote Whitcomb and Morris (1961) vage 213. 


. . . After all, any real knowledge of origins 

or of earth history anticedent to human histor- 
ical records can only be obtained through divine 
revelation. Since historical geology [or stellar 
origin], unlike other sciences, cannot deal with 
currently observable and reproductible events, 

it is manifestly impossible ever really to prove, 
by the scientific methods, any hypothesis relat- 
ing to prehuman history. 


This idea applies directly to the problem of binarv star 
evolution and origin. Experiments cannot be set up and 
conducted that will prove or disprove any case of binary 
origin. Only models of specific sets of conditions can 
be made and compared with observed binary stems. The 

"best" model which can be found has the fewest number of 


ad-hoc assumptions and predicts the greatest number of 


observed phenomena. 


CHAPTER IT 


PRE-MAIN- SEQUENCE BINARY MODEL 


The masses and periods for eleven observed close 
binary systems were selected from data tabulated by Gian- 
none and Giannuzzi (1974). The stellar structure is rep- 
resented by a polytrope of index n=3 for the radiative 
solution. A polytcope of index n=1.5 is vsed fer the Hayashi 
solution. The Hayashi solution is represented by the line 
A-B on the Hertzsprung-Russell diagram shown in Figure l. 
The radiative solution is line B-C. Both Soluvions ar¢ 
calculated for each model. The Hayashi solution is selected 
as long as it is above the radiative solution. After the 
solution is determined, the time step is calculated by an 
iterative procedure described in a later section. The mod- 
el calculations are stopped when the primary reached the 
main-sequence. 

The variables and symbols used throughout this chapter 


are listed below: 


T Kelvin-Holmholtz contraction time (years) 
At Time step between motels 

AR Radial contraction between models 

hk Stellar radius 

M Stellar mass 

ie Radial distance from Stellar Cenver 

fe) Density 

Pp Pressure 

oe Temperature 

a Radiation constant 
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FIGUPE 1. Schematic diagram of the evolutionary trecks 
on the Hertzsprung~Russell diagram for a given 
Mass. Hayashi solution line # 8). the radia= 


A 


tive solution line B - C. 


Cc Velocity of light 

Tag¢ ELtective temperature 

Y Ratio of specific heats 
n Polytropic index 

é Polytropic radius 


El Polytropic radius at surface of star 
K Constant in polytropic relationship 
K Opacity 

fo} Stefan-Boltzman Constant 

H Mass of hydrogen atom 

L Total Stellar Lumincsity 

Ln Luminosity due to nuclear sources 
lig Luminosity due to gravitational contraction 
Ly Luninosity due to the Hayashi Track 


A Separation between the centers of the binary 
components 

G Gravitational Constant 

& Total energy generéiion rate 

En Energy generation rate fer nuciear sources 

Eg Energy generation rate for gravitational contrac- 
tion 

R Gas constant 


A. Kelvin-Helmholtz Time Scale 


The Kelvin-Helmholtz contraction is known as the 
total pre-main-sequence phase of stellar evolution. The 
source of energy is derived totally from the gravitational 
potential. The time for this phase can be estimated by 
the following relationship: 


t= 1g (M/Mg) 
i/u 


fe) 
Where 19 = 2.75 X 107 years and M, is the solar mass, L 
is the solar luminosity. 

A fraction of this contrection time was used as the 


first. estimate for time stes in “he model calculations. 


The first estimated time step was found from the following: 


At = t (AR) 3.17 x 10’? seconds. 
R 
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B. Pre-Main-Seauence Mass Transfer 


The study of pre-main-sequence mass transfer shows 
that a binary will evolve rapidly from an over-contact 
configuration to a stage where both stars will fill their 
Roche Lobes. Mass transfer during condensation has been 
proposed, e. g. Wood (1962). He described a Formatian 
mechanism for a close binary in which nucici are poor in 
hydrogen compared to the mixture from which they are con-~ 
densing. Whenever two condensing nuclei of comparable 
size are present, the formation of a double star, rather 
than a single star or planets may occur. The difference 
between the formation of a double star and a single star 
is that the double star is constrained by the inner critical 
equipotential surface known as the Roche Lobe. Wood argues 
that the more massive component would continue to attract 
the most gas. 

Hayashi and Nakano (1965) and Nakano, Okyama and 
Hayashi (1968) proposed a method for the formation of a 


single star that could lead to the formation of a close 


binary. They proposed that the free-fall collapse of an 
opaque protostar is stopped by the bounce 3f a central 
core. For a protostar of one solar mass chis takes ap- 


proximately ten years. The bounce causes a shock wave. 


ae 
The subsequent propagation of the shock wave toward the 
stellar surface causes a sudden flare-up. The star finally 
settles down to a quasihydrostatic equilibrium state and 
rapidly moves toward the Hayashi Track. If the protostar 
condenses with sufficient rotational angular momentum, the 
above theory will predict the formation of a rapidly rotat- 
ing core. When the bounce shock occurs it could be non- 
spherical. The result may be the formation of a pre-main- 
sequence binary system by fission. 

Yamasaki (1971) studied the case of a binary system 
immediately after fission of a collapsing protostar. He 
suggested that a substantial fraction of the matter would 
be crowded out of the critical lobes and would drift 
around the components. The total orbital angular momentum 
of the confined mass must be above the critical value to 
cause rotational instability. He assumed that the balance 
of the mass outside of the Roche Lobes would be lost from 
the outer Lagrangean point. The loss rate is a parameter 
of the problem. If only one component overfilled its Roche 
Lobe, then the whole excess mass was assumed to fall on 
the other component. It does not escape from the system 
provided that the other does not fill its Roche Lobe. The 
model calculations were based on a polytropic star with 
n= 1.5. This represents fully convective stars contract- 
ing on the Hayashi Track. The initial mass for a 
culations was held constant at 13.5 solar masses. In all 


cases of stable mass transfer rates and mass loss rates, 


12 

the systems detached within several hundred years. For a 
critical value of mass exchange rate the system would not 
recover from the instability and form a single star. 
The final configurations of these models still had separa- 
tions from 2 to 4 A. U. Although these separations are 
a minimum, they do not represent observable close binary 
systems. 

The writer evaluated the case of pre-main-sequence 
mass transfer from the secondary to the primary. Table l 
shows the initial and final configuration after mass 


transfer. 


TABLE 1 


INITIAL AND FINAL CONFIGURATION OF MASS TRANSFER MODEL 


INITIAL FINAL 
M1/My 3.0 3.128 
M2/M, 5 R372 
Period (days) 3.0 3.459 
Separation (A. U.) 0.067 0.074 


The mass transfer rate is chosen to be 1 X 107? solar 

mass per year. This is representative of observed rates 

in semidetached binary systems, (e. g., Hall and Neff; 
1976). Table 2 shows the calculated models of the secondary 
during mass transfer. Column 1 is the model number. Column 
2 is the age of the model in seconds. Column 3 is the 
logarithm of the luminosity in solar units. Column 4 is 

the logarithm of the effective temperature. Column 5 is 


the ratio of the radius to the Roche Lobe. Finally, column 


6 is the radius in solar units. The mass transfer stops 
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when the secondary is smaller than its Roche Lobe. The 
initial radius of the secondary is arbitrarily chosen to 
be 20% greater than the Roche Lobe. Conservation of 
angular momentum and of mass requires the binary system 
to increase in separation during mass exchange. Therefore, 
the Roche Lobes of each star increase as the mass ratio 
increases. The mass transfer stops within 28 years for 
this case. This is similar to the results found by 
Yamaski (1971). The initial configuration chosen by the 
writer is immediately after the binary components form a 
minimum contact binary. Both components will fill their 


Roche Lobes initially. 


TABLE 2 


MODELS OF THE SECONDARY DURING MASS TRANSFER 


MODEL TIME LOG(L/LSUN) LOG(Tprr) ROCHE ER R/RSUN 
1 2.506E+07 4.8421 3.2855 1419 5.50 
3 6.452E+07 4.7893 3.2888 pe 5.45 
5 1.044E+08 4.7364 3.2920 1.16 5.39 
7 1.446E+08 4.6834 3.2953 Lct4 5.34 
9 1.855E+08 4.6304 3.2985 1212 5.29 

Li 2.310E+08 4.5772 3.3018 belt 5.23 
13 2.830E+08 4.5239 3.3051 1.09 5.18 
15 3.423E+08 4.4705 3.3083 1.08 5.13 
a9 4.101E+08 4.4170 3.3116 1.06 5.08 
19 4.876E+08 4.3634 3.3149 1.04 5.03 
21 5.761E+08 4.3097 3.3181 1.03 4.98 
23 6.773E+08 4.2558 3.3214 1.01 4.93 
25 7.930E+08 4.2019 3.3247 1.00 4.88 


27 9.254E+08 4.1478 3.3279 233 4.33 
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C. Basic Assumptions for the Binary Model 


The proposed model represents a pre-main-sequence 
binary star system. The initial configuration introduced 
is constrained by the Roche Lobes. The basic assumptions 
which must be made are as follows: 

1. The first law of thermodynamics requires mass 


and energy to be conserved both for the in- 
dividual stars and the binary system as a whole. 


2. The second law of thermodynamics reguires that 
the entropy for any closed system be greater 
than or equal to zero. This means that there 
is no external source (or sink) of energy 
available to do work. The system dynamics 
are determined from the initial conditions. 


3. The binary system obeys conservation of angu- 
lar momentum. Circular, Keplarian orbits were 


assumed. 


4. The stars are spherically symmetric. Thus, we 
are assuming that the stars are not affected 
by rotation during the model calculations. 


5. The stars are in quasihydrostatic equilibrium. 
This assumption is required to calculate each 
model. A contracting star is not, by defini- 
tion, in hydrostatic equilibrium. Yamasaki 
(1971) showed each contraction occurs over a 
time step which is long compared to the dyna- 
mical time scale of readjustment. Thus, the 
approximation of hydrostatic equilibrium is 
good. 


D. 


Equations of Stellar Structure 


LS 


four basic equations of stellar structure: 


ff 


Where M(r) and L(r) are the mass 


distance r. 


hvdrostatic equilibriun 


continuity of mass 


therma | 


for radiative 


equilibrium 


for convective 


equilibrium 


equilibrium 


and luminosity at radial. 


Along with these four differential equations we must 


have the following relations: 


P 


(a) 


P(p, 


T, 


T, 


chemical 


chenical 


chemical 


composition) 
zomposi on) 
composi yn) 


the equa- 
tion of 
state 


the energy 
ceneration 
rate ver 
unit mass 
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The boundary conditions for this problem are as fol- 


lows: 
At r = O, M(r) = O 
L(r) = 0 
At r =R, T = T.¢e 
P= 0 
M(r) =M 


A polytropic star is one which obeys an equation of 


the following form: 


Since the pressure is an explicit function of density only, 
polytropic stars are determined by the equation of hydro~ 
static equilibrium and the equation of continuity of mass. 
3y combining these two equations and using dimensionless 
variables, one differential equation defines the structure 


of a polytropic star: 


1_d_ (,2 a0)_ _ 9h 
Ee ae dé 
This is known as the Lane-Emden equation. The vari- 


ables are defined as follows with the subscript c indicat- 


ing the core valves. 


a ce 
Po 
6 =f 
Te 
Oe ae, 
r 
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Where ry, is known as the Emden unit of lencth: 


L 
ry, = [ie Poy? 


4mGp,2 


Le 


In addition, we have 


= Ey 


The polytropvic model vrovides a fair approximation to the 


structure of certain types of real stars. This idealized 


model is often useful in qualitative and in rough gquantitat: 


discussions and aid considerably in gaining an overall] in- 


sight into the structure of sStars,;(e. g. Cox and Guili; 1968} 


Hayashi (1961) proposed a vhase of gravitational con- 
traction of a convective star known as the Havashi Track. 
This is shown in Figure 1 as the nearly vertical line, A-3, 
on the Hertzsprung-Russell diagram. The convection may 
result from the large opacities caused by ionizing hydro- 
gen. The adiabatic relation between the pressure and 
temperature 1S a good approximation for most of the mass of 


completely convective stars. 


Pp = xp2 e 5 


Here, the ratio of svecific heats, y, iS constant and equal 
to 5/3. Even including the pvhotosyshere, this adiabatic 


relationship i183 a good approximation. 


1) 
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For stars in convective equilibrium, the adiabatic 

temperature-pressure relation can be rewritten in terms 
of Schwarzschild dimensionless variables as 


2-45 225 


“e) 
J 


{(4u/G)R4 M2} P and t = {(r/G)R/uM}T 


= 

=a 
0) 

K 

0) 

ue) 
H 


and E 4nK(u/R) 2°? qi-? ml +> Ri-S 


If the convective envelope extends all the way to the cen- 
ter of the star, then a limiting value of ™ for completely 


convective stars can be found to be: 
E = 45.48 


Hayashi described the "Forbidden Region" for a star 
of a given mass to be to the right of the line A-B. [ine 
A-B represents the surface solution for the luminosity 
based on the limiting value of E. A star finding itself 
in this forbidden region would not be in hydrostatic 
equilibrium and would have to adjust itself to another 
structure on roughly a "free-fall" time scale. 

The luminosity for fully convective gravitationally 
contracting stars, can be found by solving for radiative 
transfer in the photosphere. Cox and Giuli (1968) found 


a solution for the luminosity for stars with 


L>ly 


to have a strong dependence on the effective temperature in 


the following form: 
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Bl, B2,,B3,, B4 
Lac u M Toff 


By assuming complete ionization of hydrogen at the base of 
the photosphere, the equation for the luminosity has been 


shown to be: 


La C790 - It. 56 aoe 


Schatzman (1963) found the solution for the luminosity in 


= 


terms of solar type population TI chemical composition and 


mass to be: 


_1e 2 
Myo? ie a (Equation 1) 


Where M, is the stellar mass in solar units. 
F, The Gravitational Contraction 


For a pre-main sequence star, the only significant 
sources of energy come from the gravitational potential. 
The relationship between the gravitational contraction and 
luminosity can be determined from the equation of thermal 


equilibrium: 


L 
= Anepr2 


By using the ideal gas law, allowing the ratio of specific 


heats, y, to be equal to 5/3, and defining €= £n + €g, the 


equation for thermal equilibrium becomes 


20 
aL 


= Anqr2of aS 2IF. a. ae j 
dr DLEn (= 
2° dt 573 


By uSing the equation for hydrostatic equilibrium the fol- 


lowing relationship holds: 


32/34, Py 2.32 La 
Ze Ge 35/8 20 R dt 
Substituting, €, = 0 into the equation of thermal equilib- 


rium we get the following: 


Now this relation is rearranged in the form of polytropic 


variables: 


+1 1 dR 


Finally, the a being rewritten in different form and by 
noting that RM = €] Ary, the relation becomes 


a = - 6TrpP, €2entl AR 


Ey A 


Where AR is the radial contraction and At is the time step 
for the contraction, the total contribution to luminosity 
From a gravitational contraction of a polytvope is Found by 


integration: 


ZL 


1 
Lie 6 TP orn? AYn fs r2 pnt dé (Equation 2) 


where a = 3.69375 f£0r n = 1.5 and = 6.89685 for n = 3.0. 


This relation for the luminosity has the term Arn, The 

Ar, represents the polytropic radial sections ae At 
is the corresponding time step. Ar, 128 selected simply to 
optimize computing time. The time step, on the other hand, 
must be determined from the Virial Theorem and conservation 
of energy. 

The Virial Theorem simply states that half the energy 
made available by a contraction causes an increase in the 
stellar temperature, while the other half is radiated away. 
The potential energy Q is given by: 

R 


Q= + =} Sf Vidi. 
re) 


For a polytrope the potential energy is given by: 


o- = (25 ar 
: 5-H R 


If only half the potential energy is available for the 


luminosity, then the time step, At, is determined by 


where <L> 1S the average value of tne luminosity for corres- 
g 


ponding time step. Therefore, we have 
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22s 
itl Ry 


(Equation 3) 


and <L> = (Ly4 7 + L3) 


where the index i corresponds to the ith model. 

The problem then is to determine, by an iterative 
procedure, selfconsistent values for the new value of the 
luminosity Lj4, and the time step tj,) - ty. An initial 
guess for the time step was obtained from the Kelvin- 
Helmholtz time scale for contraction onto the main-sequence. 
The estimate was obtained by proportioning the time scale 
for the radial contraction step. An estimate for the lumin- 
osity, L', is calculated for the polytrope. An estimate 
for the time step, At', is calculated from L'. Then using 
At' a new value of the luminosity, L", is calculated and 
compared to L'. If L' is greater than L", then a new in- 
creased time step, At", is estimated. Likewise, if L' is 
less than L", then the new estimate for the time step is de- 
creased. This can be written as follows: 

At = At" +/At' - at"|” 
where -"- + "was used 1f L"<L" or " - was used if L" oL'. 
The primed values are the old estimates and the double prime 


indicates the new estimate. After L" and L' approach each 


Z3 


other to the desired accuracy then L" is set equal to Lj;4}.- 
The time step, At", is stored for that model as At and be- 
comes the first estimate for the time step for the next 


model. 


G. Nuclear Energy Sources 


If the assumption is that no contraction takes place 


then the equation for thermal equilibrium becomes 


By using the appropriate polytropic variables and the nuclear 
energy generation rate fn for the proton=proton reaction; 


the luminosity is 


E 
4 2 2q2nt4 é 
Ly = 47K, 0," T, EX rere” a& (Equation 4) 


where €4 = 10 x 10739 ergs/sec/gm and X is the hydrogen 
fraction of stellar material. 

The nuclear luminosity is calculated and compared to 
the luminosity of the "zero-age" main-sequence star of 
corresponding mass. The model calculations ended when Ly 
approached the main-sequence value. 

The time that a particular component of the hinary 
remained on the main-sequence is considered to evaluate 
the systems. Generally the primary should reach the main- 


sequence before the secondary. It is of interest whether 
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the primary has enough time to evolve off of the main-se- 
quence before the secondary can reach the main-Sequence. 

A logarithmic regression is »verformed on data published by 
Iben (1965), and Strothers (1963) (1964) (1966). The fol- 
lowing Table 3 represents mean published values for the time 


on the main-sequence tug, Novotny (1973). 
TABLE 3 
TIME ON THE MATN-SEQUENCE 


M/My Log M/M, log tys 
60 1.7782 6.4914 
45 1.6532 6.5563 
30 1.477] 6.6721 
LS 1.1761 7.0043 
9 -9542 7.3242 
5 -6990 7.8088 
3 4771 8.3443 
2.25 23522 8.7259 
1.5 -L761 9.2966 
1.25 .0969 9.6053 
1.0 0 10.0086 
0.5 -.3010 10.4771 
0.25 -.6021 10.8541 


The following relationship can be used to represent the 
above logarithmic values: (The maximum error is less 


than 5% over the entire range of masses.) 


log tya = (-2.00998) log (F) + 9.62409 


ae) 


MS 


In all cases evaluated, the secondary has sufficient time 
to contract onto the main-seguence before the primary 


evolves away. 
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H. Initial Conditions 


With the masses of each comoonent and their veriod 


being free parameters, the Roche Lobe for each star was 


determined by the program. The following relation, given 


by Paczynski (1971) for the Roche radius, R,, was used: 
Rr = A (0.38 + 0.2 log (My /Mg) ) 


The Roche radius was determined for cach star and was 
taken to be the initial stellar radius for the First model 
of each star. The program then calculated the sequence of 
contraction models. 

The chemical composition is chosen to be of popula- 
tion I type and assumed to be the same for each star. The 
hydrogen fraction, X, is selected to be 0.750, the helium 
fraction, Y, selected to be 0.224, and the heavy metal 
fraction, Z, is selected to be 0.026. This selection is 
arbitrary but is chosen to allow comparison with the work 
of Iben (1965). 


The total mass for each model is between 2.5 and 6 


solar masses. Eleven close binary systems are chosen from 
Glannone and Giannuzzi (1974). The parameters are listed 
in Table 4. Column 1 is the name of the System. Column 2 


is the period in days . Column 3 is the mass in solar units 
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for both the primary and the secondary. Column 4 is the 
logarithm of their luminosities in solar units. Finally, 


column 5 is the logarithm of their effective temperatures. 


TABLE 4 


OBSERVED SEMIDETACHED BINARY SYSTEMS 


SYSTEM PERIOD M/My LOG (L/Lo) LOG (TprrF) 
TV Gas 1.243 3.10 2.070 4.029 
1.39 1.028 3.787 
IN Aur 1.245 2.97 1.999 4.029 
0.89 0.203 3.679 
U Ceo 2.493 3.19 2.066 4.079 
1.53 9.875 3.678 
AI Dra 1.199 2.18 1.467 3.982 
1.03 0.659 3.756 
6 Lib 2.327 2.96 1.976 3.982 
1.31 0.779 3.675 
8 Per 2.867 3.15 2.082 4.079 
0.74 0.797 3.696 
U Sge 3.381 4.27 2.593 4.134 
1.60 0.573 3.567 
V505 Sag 1.183 2.22 1.501 3.958 
1.18 0.414 3.693 
X Tri 0.972 1.72 1.044 3.947 
1.00 0.000 3.642 
TX U Ma 3.063 3.13 1.997 4.079 
0.99 1.088 3.746 
WU Mi 1.701 2.68 1.824 3.947 
1.19 0.541 3.681 


I. Description of the Model Program 


The program calculates both stellar model sequences 
Simultaneously. This allows comparison of time steps and 
Structure of each component when vossible. The initial 
mass of each component and the period of the binary are 
barameters of each model sequence. Each observable system 


can be evaluated in turn. 
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Figure 2 represents a flow chart for the pre-maimse- 
quence model program. By using the Kelvin-Helmholtz con- 
traction time, an estimate for the time stev, At, is made 
for the initial radial contraction AR. The initial value 
of the luminosity, L', is arbitrarily chosen from the mass- 


luminosity relation for main-sequence stars: 


log 42.) = 44, leg (4 = 001 
i Mo 


The radial contraction, Arn, wes chosen in alli cases to be 


(Arn) 544 = 9.005 Ri 
El 
A polytropic model is calculated for both n = 1.5 and n= 3. 
By uSing Equation 2, the first value for Lg is calculated. 
The effective temperature is calculated by using the follow- 
ing relation: 


4 


= L 
T os 
eff Ato R2 


By using Equation 1, the Hayashi Track luminosity, Ly, is 
found and compared to the radiative luminosity, Lg. If 


Lg<ly, then L, is chosen; if not;then Ly is chosen. Now 


g 
the time step can be calculated from the Viral Theorem 
using Fgquation 3. If the old value for the luminosity 


L 18 equal to Lb” within suriicient.accuracy, them the 


procedure is stopped; if not, the time step is adjusted 


rs 
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FIGURE 2. Flow chart of the program for the pre-main-sequence 
model 
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as described previously and the procedure starts over 
again. After the error between L" and L' is less than 1 
part in 10°, the values for R, Tor¢, L, and At are stored 
for that model. If the program limits have been met, L, 
is large enough; the program stops; if not a radial con- 
traction is made and the procedure is restarted. 

Appendix A lists the main program, all subroutines 
and polytropic constants. The language is Honeywell's 
Mark III version of General Electric Company's BASIC. 

A test model is compared with the evolutionary tracks 
published by Iben (1965). The primary for the test model 
is 3 solar masses, the secondary is 1.5 solar masses, ané 
the period is 3 days. Initially both stars are filling 
their Roche Lobes and are on the Hayashi Track. Figure 3 
shows the plots of the test model and Iben's pre-main- 
sequence tracks for single stars on the Hertzsprung-Russell 
diagram. The initial point for Iben's models is labeled 
"QO". The points 1] through 5 can be used for age compari- 
sons with the test model. Table 5 shows the ages at these 
points labeled on the evolutionary tracks in Figure 3 for 
each star. The same numbering is used for the test model. 
The test model tabulation for both the primary and the 
secondary can be found in Appendix B. Each table has 6 
columns. Column 1 is the model number after the initial 
model. Column 2 is the accumulated time in seconds after 


the initial model. 
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FIGURE 3. Evolutionary tracks of a pre-main-sequence test 
model for a binary (@= 3M5, @ = 1.5 Mo) and the 
evolutionary tracks for single stars published 
by Iben (1965). 
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TABLE 5 


COMPARISON OF AGE (106 yEARS) FOR THE TEST MODEL 


PRIMARY = 3 M SECONDARY = 1.5 “% 
AGE POINT IBEN TEST MODEL IBEN TEST MODEL 
1 0.034 = 0.23 0.205 
Za 0.208 0.06 2.36 Peso 
3 0.763 0.747 5680 3.483 
4 6 Rn ibe 1.202 7238 4.962 
5 1.250 bees eo al 8.62 5.949 
6 1.465 1.432 10.43 Saas 
7 1.741 SS a A ee nS 
8 2.514 as Es.21 soot 


Column 3 is the logarithm of the luminosity in solar units. 
Column 4 is the logarithm of the effective temperature. 
Column 5 is the ratio of the model radius to that of its 
Roche Lobe. Finally, column 6 is the radius of the model 
in solar units. 

These ages and evolutionary tracks of the test model 
show reasonable agreement with Iben's results. Up to 
point 5, the age of a component of the test model compares 
with the age of Iben's models. Suggestions for future im- 
provements and changes are discussed in Chapter IV. The 
major differences in ages occur near the main-sequence. 
These differences are caused by the limitations of a voly- 
tropic model to adjust for the core's changing chemical 
composition. In particular, Iben followed ‘he depletion 
of Cl2 in detail. Normally, cl2 is mainteined in equilib- 
rium for the CN reaction chain. However, during the pre- 


main-sequence contraction the CN reaction has not started. 
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Only very near the main-sequence is there sufficient ther- 
mal energy in the core to initiate this reaction. cl2 
must be depleted to form the equilibrium concentrations 
of the "secondary" elements in the CN reaction chain. 
After equilibrium is reached, the concentration of cl2 
remains constant. This resulted in slowing down the con- 
traction near the main-sequence as Cl2 wag depleted in 
the core. The test model is polytropic and cannot follow 
detail changes in stellar structure. It is intended as 
a survey tool to determine whether or not a binary system 


is pre-main-sequence. 


CHAPTER IIi 


RESULTS OF THE NUMERICAL CALCULATIONS 


The age computed by the writer's model is the time 
after both components filled their Roche Lobes. It may 
not be the same as would be calculated for a single star 
contracting onto the main-sequence. 

Appendix C lists the numerical values for each of 
the pre-main-sequence models calculated as described in 
the previous chapter. The Appendix C contains 22 tables. 
Each table is identified by the name of the binary system. 
Each table has 6 columns. Column 1 is the model number 
after the initial model. Column 2 is the accumulated time 
in seconds after the initial model. Column 3 is the 
logarithm of the luminosity in solar units. Column 4 is 
the logarithm of the effective temperature. Column 5 is 
the ratio of the model radius to that of its Roche Lobe. 
Finally, column 6 is the radius of the model in solar units. 

The computed binary model and an observed binary sys- 
tem are compared by age and position on the Hertzsprung- 
Russell diagram. The ages of each component ‘aust match and 


the observed binary must fall near the evolutionary tracks 
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of the binary model for the systems to be pre-main-sequence. 


Individual systems are discussed in the following sections. 
A. TV Cassioneia 


It appears unlikely that TV Cassiopeia is a pre-main- 
sequence binary. The secondary is either éoee—Tuntheus ox 
it has too great an effective temperature. This is in con- 
tradiction to the results published by Field (1969) for 
TV Cassiopeia. Field determined the pre-main-sequence 
nature of the system by estimating the radii and age of 
each star. The estimate is made by interpolation of pub- 
lished models of a single star. 

Figure 4 shows the plotted contraction of both com- 
ponents of a pre-main-sequence model for TV Cassiopeia on 
the Hertzsprung-Russell diagram. The observed luminosity 
and effective temperature for both components are plotted 
also. The problem of the over-luminosity of the secondary 
is clearly shown. 

If the age of the primary can be determined from the 
model, it is approximately 6.3 XK 10° years since the in- 
itial model. To estimate the age of the secondary, either 
the radiative solution or the Hayashi solution must assume 
to be correct. Comparing the luminosity with that of the 
model secondary on the Hayashi Track yields an age of 1.5 
x 10° years. If the radiative solution is assumed, the 
age is 3.1 ¥ 105 years. ‘his pre-main-sequence model does 


not fit the observed system for either case. 
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FIGURE 4. TV CASSIOPEIA (A = 3.1 Mo, A = 1.39 Mo). The 
evolutionary tracks of the pre-main-sequence 
binary model is represented by the circles 
(primary @e, secondary 0). 
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B. IM Aurigae 


Figure 5 shows the plotted contraction of a pre- 
IMain-sequence model for IM Aurigae on the Hertzsprung- 
Russell diagram. The observed stellar luminosity and 
effective temperature for both components also are plotted. 
The observed system falls close to the contraction paths 
of the pre-main-sequence model. However, the "ages" do 
not match. 

The age of the primary seems to be 3.1 X 10> years. 

The age of the secondary appears to be between 9 X 105 


and 3.1 X 10® years according tc this mo-2l. 
C. U Cepheil 


Figure 6 shows the plotted contraction of a pre-main- 
sequence model of U Cephei on the Hertzsprung-Russell dia- 
gram. The observed luminosity and effective temperature 
for both components also are plotted. This system appears 
to be pre-main-sequence. 

First of all, the observed system Falls close to the 
contraction path of the model. This is true for both com- 
ponents. Secondly, the ages match within 20% at 3.15 X 


10° years if the secondary is on the Hayashi Track and has 


2) 


a slightly different chemical composition. TIben (1965) 
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FIGURES 5. TM AURIGAE (A = 2.97 Mo) A = 0.69 Mo) : The 
evolutionary tracks of the pre-main-seguence 
binary model is represented by the circles 
(primary @, secondary 0). 
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FIGURE 6. U CEPHET (a = 3-19 My, 4 = 1-53 M,). The evolu- 
tionary tracks of the pre-main-sequence binary 
model is represented by the circles (primary e, 
secondary 0). 
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showed a slight decrease in the metals composition would 
move the Hayashi Track toward higher effective temperatures. 
This would match ages and positions on the Hertzsprung- 
Russell diagram. This model predicts that U Cephei is a 
pre-main-sequence binary star system with population type 


I chemical composition for both components. 


D. AI Draconis 


Figure 7 shows both the path of the pre~-main-sequence 
contracting model and the observed binary system AI 
Draconis plotted on the Hertzsprung-Russell diagram. The 
system does not appear to be a pre~main-sequence object 
according to this model. This problem seems to be the 
same as TV Cassiopeia. The ages do not match. The secon- 
dary component of the observable system is over-luminous; 
or the effective temperature is too high for it to be con- 


sidered a pre-main-sequence star for this particular mass. 
Be. 6 Libre 


Figure 8 shows the pre-main-sequence contraction model 
for 6 Librae plotted on the Hertzsprung-Russell diagram. 
The observed binary system may possibly be pre-main~ 


sequence if it is assumed that the two stars have aitferent 


chemical compositions. This would follow the description 
by Wood (1962) of the formation of a binary with different 


compositions. As stated before, Ibsen (1965) has shown that 


such a change in composition could increase the effective 
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tionary tracks of the pre-main~-sequence binary 
model is represented by the circles (primary @, 
secondary o). 
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temperature of a star on the Hayashi Track. With this in 
mind, the ages of both components are in agreement. The 
age after both components filled their Roche Lobes is 


approximately 8.2 X 105 years. 


F. 8 Persei 


Figure 9 shows the observed binary 8 Persei plotted 
along with the path of a pre-main-sequence contracting 


model. This system does not fit the model for a pre-main- 


sequence binary star system. It suffers from the same 
problems as TV CaSsiopeia and AT Draconis. ‘The secondary 
is not represented by this pre-main-sequence model. The 


age is also a problem. If the chemical compositions are 
different enough to fit the model, the age for the secondary 
would be younger than the primary by a large factor. This 
is caused by the decrease in contraction time as the lumin- 
osity increases. The luminosity must be much higher for a 
given effective temperature for the secondary to fall on 


the Hayashi Track. 


G. U Sagittae 


Figure 10 shows the plot of the binary system U 


Sagittae on the Hertzsprung-Russell diagram. Comparing 
the evolutionary tracks of the model with tne observed 
system, the binary could be pre-main-sequence. The ages 


of both stars appear to match at approximately 3.78 X 10° 
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FIGURE 10. U SAGITTAB (A = 4.27 Mo, A = 1.60 Mo). The 
evolutionary tracks of the pre-main-sequence 
binary model is represented by the circles 
(primary @, secondary 0). 
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years. An increase in metals content of the secondary 
would move the Hayashi Track toward lower effective tem- 
perature. U Sagittae is probably a pre-main-sequence ob- 
ject. Therefore, another alternative is that U Sagittae 
may not have a fully convective pre-main-sequence secondary. 

The secondary falls in the forbidden ones described 
by Hayashi. If this is actually the case, then the star 
is not in hydrostatic equilibrium, and its structure would 
be changing rapidly. Larson (1969) described pre-main- 
sequence contractions that may not be on the Hayashi Track. 
These evolutionary tracks had rapidly varying luminosities 
as a function of effective temperature, caused py different 


depths of the convective envelope. 
Hs.  VS05 Sagittarii 


Figure 11 shows the paths of a contracting pre-main- 
sequence model plotted on the Hertzsprung-Russell diagram. 
Both components of the observable binary system V505 Sagit~ 
tarii are plotted along with the model. By using the same 
argument that the components could have slightly different 
chemical compositions, the system could be pre-main-sequence. 
The age of the system would be approximately 1.89 X 106 


years. 
I. X Trianguli 


Figure 12 shows the paths of the contracting binary 


model and both components of the binary X Trianguli plotted 
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V505 SAGITTARII (A = 2.22 M,, A= 1.18 M). 
The evolutionary tracks of the pre-main- 
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sequence binary model is represented by the . 


circles (primary », secondary 0). 
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FIGURE 12. 
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X TRIANGUILI (4 = 1.72 M,, A= 1.00 MQ). The 
evoltuionary tracks of the pre-main-sequence 
Dinary model is represented by the circles 

(primary @, secondary o). 
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on the Hertzsprung-Russell diagram. The observed binary 
system falls close to the pre-main-sequence path of the 
model. The ages, however, do not quite match. The dif- 
ference in the ages is a little less than a factor of two. 
This could possibly be caused by a difference in chemical 
composition. If we allow for a slight increase in effec- 
tive temperature, caused by a decrease in heavy metals, 
then the time of contraction would decrease. This is re- 
quired by conservation of energy and the Virial Vheoren. 
With these assumptions the binary X Trianguli could pos- 


sibly be a pre~main-sequence object. 
J. TX Ursae Majoris 


Figure 13 shows the observed binary TX Ursae Majoris 
plotted on the Hertzsprung-Russell diagram. The contract- 
ing pre-main-sequence model is plotted also. The system 
does not appear to fit the model. The age would have to 
be =9.5 X 109 years as determined by the primary. This 
yields the same problems that were found with 8 Persei and 
others. The secondary is over-luminous for a pre-main- 


sequence star with this particular mass. 
K. W Ursae Minoris 


Figure 14 shows the observed contact binary system 
W Ursae Minoris plotted on the Hertzsprung-Russell diagram. 


The path of the contracting pre-main-sequence model also 
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FIGURE 14. W URSAE MINORIS (A = 2.68 M4, A = 1.19 My). 


The evolutionary tracks of the pre-main- 
sequence binary model is represented by 
the circles (primary e, secondary 0). 
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is plotted. Both components of the observable system are 


relatively close to the path of the model. The exciting 
part of this analysis is that the ages for the pre-main- 
sequence model match at about 1.87 x 106 years for a sys- 
tem slightly different from the observed system. Iben 
(1965) showed that an increase in metals could account for 
the increased luminosity of the primary on the radiative 
path. He also showed that a decrease in metals could in- 
crease the effective temperature for the secondary on the 
Hayashi Track. The results of this model cannot be dis- 


counted concerning the possible pre-main-sequence nature 


of W Ursae Minoris. 


CHAPTER IV 


CONCLUSION AND DISCUSSION 


This binary model allows discovery of pre-main- 
sequence systems. Initially both stars must be forced 
into their inner equipotential surfaces rapidly. Evolu- 
tionary tracks calculated from this initial configuration 
define a new means for evaluating actual binary systems. 

The major conclusions cf this evaluation are as 
follows: 

l. Pre-main-sequence close binary star systems can 
be found by comparing observed binaries with this 
model. Some observed closed binary systems are 
found to fit this pre-main-sequence model. 

2. The chemical composition of the two components 
of the pre-main sequence binaries appeared to be 
slightly different. Even though the actual dif- 
ferences in composition were not quantified, 
both stars could still be considered population 
type I. 

3. There does not seem to be a correlation between 
the orbital angular momentum, ages, or mass ratios 
for the pre-main~sequence binary systems. 
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Table 6 summarizes the evaluation of the eleven 
binaries. The distribution falls fairly evenly into three 
groups, those binaries that are likely to be pre-main- 
sequence, those that are possibly, and those that are un- 
likely to be. Improvements in the model could allow better 
determination of the pre-main-sequence nature of a binary 
system. A more realistic model, instead of assuming a 
polytropic structure would guarantee more accurate struc- 
ture and age determination. It would allow for changing 
the chemical composition of both stars. The close binary 
systems under consideration most likely are in synchronous 
rotation. The periods are approximately a few days. This 
rotation has been neglected in the calculations, but it 
surely would slow the gravitational contraction to some 
extent. 

Roxburgh (1967) evaluated the effects of rotation and 
magnetic fields on pre-main-sequence evolution of binary 
systems. Rotation of an individual star could be uniform 
or differential. Uniform rotation requires a large vis- 
cosity of the stellar material to couple the core to the 
envelope. If the star is in convective equilibrium or has 
a sufficiently large magnetic field (>1 Gauss), then the 
viscosity of the stellar material should he great enough 
to allow uniform rotation. In this case, Roxburgh predicts 
loss of material at the stellar equator. Thus, rotating 
pre-main-sequence stars in convective equilibrium will fol- 


low a more vertical Hayashi Track. 
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"con- 


Stars in radiative equilibrium do not have the 
vective viscosity" mechanism. They are more likely to 
approach the critical angular momentum valve and fission 
may result. He concludes that contact binary systems could 
be formed by fission of a single pre-main-sequence star 
on the radiative track. For binary systems that are formed 
before they reach the radiative track, the rotational 
angular momentum should not reach the critical value. 

Both components should be in synchronous rotation and most 
of the angular momentum would be orbital. 

Bierman and Hall (1976) and Oliver (1978) cvaluated 
the possibility that RS Canum Venaticorum systems could be 
pre-main-sequence. Their arguments were based on space 
densities and probabilities assigned by evaluating ages. 
Because of the large number of these systems, they concluded 
that they must be post-main-sequence. This conclusion is 
valid as long as the universe is the accepted cosmological 
age. However, there is no direct observational evidence 
that these systems must be post-main-sequence. The writer 
would like to apply an appropriate pre-main-sequence model 
to RS Canum Venaticorum type systems. Bierman and Hall 


concluded that RS Canum Venaticorum system resulted from 


the fission of a post-main-seguence star. It would be in- 
teresting to model a thermal contraction of a hydrogen de- 
pleted secondary. This may represent a binary immediately 


after fission. 


56 
In studies of pre-main-sequence and post-main-sequence 
evolution, the origin of a forming system is assumed. The 
writer is assuming the origin of the binary system. 
Clayton (1978) and Larson (1978) and others have presented 
theories for stellar formation. They admit that unresolved 
problems remain. These include gas dispersion, viscosities, 
and turbulences within a collapsing gas cloud. The prob- 
lems are more severe for the formation of a binary star 
system. Many scientists have assumed the origin when mathe- 
matical modeling beyond a particular point is difficult or 
impossible. Whitcomb and Morris (1961) have proposed 
special creation beyond the known laws of physical science. 
Special creation states that God created the universe 
(Genesis 1:1). He created each object with a specific 
purpose ( I Corinthians 15:40, 41) 
. . . There are also celestial bodies and bodies 
terrestrial: but the glory of the celestial is 
one, and the glory of the terrestrial is another. 
There is one glory of the sun, and another 


glory of the moon, and another glory of the stars: 
for one star differeth from another star in glory. 
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PRE-MAIN--SEQUENCE BINARY MODEL PROGRAM LISTING 
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149 Pel TEHT 3 1.40207 290590 J47 2.04 
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MAYOET TIMc TASCE JEGERY) TAG( TONE) DACUE =O DOSE 
159 DPLTIABU HY 2 1 ,AIAD 2,90 75 eg D299 
15] 2eTARIEF] 3 1 ADAG 2,920 Al 2.92 
15? DLTILE4I2 1 4227 2.9498 744 2K 
152 2. fOSS41 32 1.4300 29294 WA Be) 
154 De PORE BS Yl oa] 2,934] 214 292 
155 2.82KE412 1 4353 U09257 J 44 2.92 
154 DL BAATHEITR 1 LAQTA WIBATA JAA 2.9% 
157 D2 ,RADEF+] 2% 1.4308 3.9207 JA5 2629 
152 a SA1ES13 1, 4406 2, 9ANT 245 Pig 
159 2,%O0R+{ 2 1.4442 3,942 2% 0 A5 2 LRA 
149 2.901 7TEH13 | KAAS 2.9429 15 2,84 
141 D2 ,AIBE47 2% 1 .44R7 29454 ~ AA D 32 
142 pe AEeIS 145500 2, 9ATO 244 gee, 
142 DL OTPE+1 2 1.4522 QALQARO ~A4 2 2an 
144 2.000F+1 3 1.4554 2.9595 ~ 44 X70 
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TABLE B-2 
PRE-MAIN-SEQUENCE BINARY TEST MODEL FOR THE SECONDARY 
MASS = 1.5 M, 


PERIOD = 3 DAYS 
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te 1, O2P9E+99 2.2592 2 27 1 23 
[a 2, NOAE4NO 2,224D 2, 37D4 2] AO} 
19 2. 2ASEF+09 2.24592 a FAR AY AAO 
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APPENDIX C 


TABLE C-l 
PRE-MAIN-SEQUENCE BINARY MODEL FOR TV CASSIOPEIA 
PRIMARY MASS = 3.10 M 


PERIOD = 1.243 DAYS 


MY NET Trte ahey Ae. TGCS sey oarNe ep Qasr 


1 1.447541] 1.598) 3, 20A7 oy ALA 
? 3, 2904R tt 1.579 2, AOA 2? ae ae 
2 A,DADE4IT 1.8} 25 2.30aN 278 ASA 
4 A.DP90ES I] 1.5144 ILANNA mel) A,54 
5 Aig 2a ey Leh] Ae SS yer 4.5) 
A O88 4 14 1.5199 2.2399 oO? 4,49 
a P.t52%S419 1.59190 2.2145 6 AAT 
R T.31eF4i2 1.5224 FoF AD -OF ALAS 
9 1.4828410 $5985 Fig Si OS AAD 
19 1.44754 12 1.5277 2, 8},94 295 ALAN 
im P.312F+I2 81.5290 2,274) OA 4,3 
12 Leet Tet? Tee i2d 46 SAO oA AYA 
13 2.14154+12 1.5342 ALADAYX 92 AL3A 
14 DLICKAE4I2 1. 52AA 3,204) 223 AWD 
15 2e*7ILEtI2 1,453R46 PQ OES 4? 4.29 
14 Ae Ihe] 2 Lea FOP 3,9292 e22 Aioet 
17 DLANNT4+12 1.5420 2 B2IN9 mes] 4 dD= 
18 2.°O5F419 1.5457 4 8325 et A232 
ho a, ) ase” Dhara 2,924] an ALDt 
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PRE-MAIN SEQUENCE BINARY MODEL FOR W URSAE MINORIS 
PRIMARY MASS = 2.68 M, 
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